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Many organisms form colonies for a transient period of time to withstand environmental pressure. 
Bacterial biofilms are a prototypical example of such behavior. Despite significant interest across 
disciplines, physical mechanisms governing the formation and dissolution of bacterial colonies are 
still poorly understood. Starting from a kinetic description of motile and interacting cells we derive 
a hydrodynamic equation for their density on a surface. We use it to describe formation of multiple 
colonies with sizes consistent with experimental data and to discuss their dissolution. 


Colony formation is a pervasive phenomenon in living 
systems and is crucial for the survival of many species 
[6]. One of the well-known examples where colony for¬ 
mation is essential are biofilms. A bacterial colony can 
grow from a single cell via multiple cell divisions mm- 
However, there is another mechanism, which relies on 
successive encounters of individual, motile bacteria, as 
also occurring in the initial stages of biofilm formation. 
This scenario of a kinetic formation of colonies domi¬ 
nates over proliferation if individuals are highly motile 
and their encounters drive the assembly of cells on a time 
scales much shorter than the characteristic cell division 
time. N. gonorrhoeae or N. meningitidis on biotic or 
abiotic substrates such as glass [7], plastic (Fig. 0^)) or 
epithelial tissue [3] are prototypical examples for such a 
scenario. Motility of these and many other bacteria origi¬ 
nates from long and thin filaments, called pili, which grow 
out the cell, attach to a substrate, retract and thereby 
actively pull the cell forward IHHII]- Pili are also used to 
mediate attractive displacements between cells laiiiHisi 
with a characteristic interaction scale given by the mean 
pili length. Colonies begin to form within thirty min¬ 
utes, which is significantly smaller than the characteris¬ 
tic cell division time-scale (A. gonorrhoeae: approx. 3 
hni). Bacterial colonies are in general reversible struc¬ 
tures. Under certain conditions, for example the lack 
of nutrients or oxygen, they can dissolve and re-colonize 
their surroundings [EHII]. Specifically, N. meningitidis 
and N. gonorrhoeae bacterial colonies have been shown 
to dissolve by effectively lowering the strength of the pili- 
mediated interaction mni]. 

However, so far, the physical mechanisms governing 
the formation and dissolution of bacterial colonies are 
poorly understood. Since motility and interactions are 
driven by active retractions of pili, fundamental concepts 
from equilibrium statistical mechanics are in general not 
applicable. The inherent non-equilibrium nature of this 
system suggests to consider a kinetic approach reminis¬ 
cent of the Boltzmann equation, which has been success¬ 
fully employed to describe the order-disorder transitions 
in several active systems far from equilibrium p!8ti28] . 

Here we propose a kinetic description as a general 
framework of how living colonies form and dissolve, which 
keeps track of the length scales and the specific proper¬ 
ties of the interactions between individuals. By a coarse- 



FIG. 1. (color online) (a) N. gonorrhoeae colonies three 
hours after sedimentation on a plastic substrate, (b) Typi¬ 
cal snapshot of state obtained from a numerical solution to 
Eq. § at large time scales. 


graining procedure we derive the corresponding hydro- 
dynamic equation and find an ordering instability for 
a choice of parameters relevant to N. gonorrhoeae. It 
belongs to a class of instabilities, where the diffusion 
constant is negative and originates from attractive pili- 
mediated interactions. As most of the parameters can be 
estimated based on available data for N. gonorrhoeae, we 
analytically compute the corresponding phase-diagram 
and the characteristic colony size that is consistent with 
experimental observations. 

Our theory can also be used to compare the effects of 
different cell-cell interactions and investigate their inter¬ 
play. We show that pili interactions are more effective 
regarding clustering than cell adhesion. Moreover, when 
both interactions keep the cells together in the colony, a 
more efficient and robust way to dissolve the colony is 
to lower the strength of pili-mediated interactions. This 
suggests that pili play an essential role not only in cell 
motility and assembly, but also in the dissolution of ma¬ 
tured colonies. Our results demonstrate that kinetic the¬ 
ory can be applied to quantify the process of colony for¬ 
mation in living systems and is able to provide insights 
about the underlying physical mechanisms. 

Kinetie Model: Our kinetic description is formulated in 
terms of the particle density /(r, t). We restrict ourselves 
to two-dimensional colonies forming on a planar sub¬ 
strate [29] , which do not give rise to swarms or swirls (see 
e.g. [22]). Therefore, the spatial coordinates r G suf- 
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fice as dynamical variables. In the absence of interac¬ 
tions cells are assumed to move across the substrate by 
pili-mediated displacements as in case of N. gonorrhoeae 
or N. meningitidis^ leading to a diffusive behaviour at 
large length and time-scales m- Interactions enter the 
kinetic description via “collision rules”. A collision rule 
IZ maps the pre-collision coordinates to the post-collision 
positions by means of the delta-functions S{'). The cor¬ 
responding kinetic equation is: 

dtf{r,t) =Cmot{r,t)+Cint{r,t), (la) 

where Cmot describes the cell motility across the substrate 

Cmot(i',i) = J Ar' [K-r /^r /(!•', i) - JCr^r' /(r,i)] (lb) 

and Cint accounts for the cell-cell interactions 

Cint{r,t) = t J dn ydr 2 W(|ri 2 |)/(ri,t)/(r 2 ,t) (Ic) 

X [5 (7^i(ri,r2) - r) + <5 (7^2 (ri,r 2 ) - r) - 2S (r 2 - r)]. 

JCr^r' denotes the transition kernel to move from r to 
r' by a retraction event of an individual pilus. We as¬ 
sume that retraction events are independent and that 
the corresponding rate is isotropic, with a characteristic 
length scale given by the pili length £pi. There is ex¬ 
perimental evidence that the pili lengths are distributed 
exponentially [30]. Therefore, we consider for the tran¬ 
sition kernel /Cr^r+b = ACo/(27r^pi) exp ( —|b|/-^pi), with 
/Co denoting the attachment rate of pili to the substrate 
and b = r' — r is the displacement resulting from an 
individual pilus retraction. 

W(|ri 2 |) characterizes the isotropic kernel for collisions 
between cells with |ri 2 | = |ri — r 2 | denoting the relative 
cell-cell distance. For pili-mediated attractive displace¬ 
ments, we consider the following collision rule: 

(ri, r 2 ) (7^l, 7^2) = (ri - ari 2 , T 2 + ari 2 ), (Id) 

where a G [0,1/2] is a measure for the strength of the 
attractive interaction. For a = 1/2, cells are maximally 
attracted and displaced to the center-of-mass coordinate 
Ri 2 = (i*i + r 2 ) /2 between the collision partners, while 
for a = 0, cells diffuse freely without interacting. Due to 
the exponential distribution of the pili lengths, the inter- 
action rate is Wpi(|ri 2 |) = 7Wo/(27r^2.) g^p (-|ri2|/4i), 
where sets the characteristic length scale for the at¬ 
tractive interaction and Wo denotes the interaction rate. 
Since pili-mediated cell-cell interactions are intrinsically 
stochastic 13 E], we introduce a non-dimensional num¬ 
ber, 7 , accounting for the number of successful binding 
and retraction events to the total number of pili-cell en¬ 
counter events. 

Coarse-graining: The isotropy of the interaction rates 
allows us to integrate Eq. Q over the center-of-mass 
coordinates R 12 leading to non-local terms (see Sup¬ 
plemental Material [31], SI). These terms are related 
to the length scales of the interactions and resemble 


a phenomenological description for the assembly of ac¬ 
tive bundles [32H34] . Since cell colonies typically ex¬ 
hibit sizes noticeably beyond the interaction length scale, 
the non-local integrands can be removed by expanding 
the particle density / with respect to the spatial coor¬ 
dinates [ism]. Truncation of this expansion amounts 
to coarse-graining beyond the interaction length scale. 
To obtain a well-defined set of hydrodynamic equations 
for the dynamics of bacterial colonies with pili-mediated 
interactions we truncate at the fourth order (see Supple¬ 
mental Material [31], S 2 ): 

dtp{r,t) = a{p)V^p{r,t) - j3i |Vp(r,i)p + k{p) VV(r,t) 
+ ^2 [V7(r,t)]'-/33 [Vp(r,i)]-V7(r,t), (2) 

where p = / • ^p- is the dimensionless density and the 
kinetic coefficients are 

a{p) = G-^ip(r,t), (3) 

^(p) = -(^2 T/^s) p(r,t), (4) 

Pi = adC2, P2 = a^d^C4/4, /^s = (ad^ + da^) C4/6; 

d = 1 — a. Note that all Pi > t). The numerical con¬ 
stants Ck are given in the Table of the Supplemental Ma¬ 
terial [31] . In Eq. we rescaled coordinates by the pili 
length ^pi, i.e. r ^ r • ^pi, leading to a rescaling of time 
t ^ tTp-/(W 07 ). We introduce the dimensionless param¬ 
eter, G = D/( 7 Wo), with D = 3/Co'^pi denoting the sin¬ 
gle cell diffusion constant. G is reminiscent of the inverse 
Peclet number and can be interpreted as a measure for 
the rate of diffusive particle transport relative to the fre¬ 
quency of interactions. In other words, given a time pe¬ 
riod between two successive collisions, G quantifies how 
much distance is traveled (on average) by diffusion with 
respect to the mean free path. 

A equation similar to Eq. but phenomenologically 
constructed appeared in the context of laminar flames 
and propagation of concentration waves referred to as 
Kuramoto-Sivashinsky equation [35] [36]. It has also been 
pointed out as an appropriate framework to study insta¬ 
bilities in growing yeast colonies [ 6 ]. However, Eq. ([^ 
is distinctively different because the kinetic coefficients 
depend on density [Eqs. © and ©]. Moreover, Eq. © 
exhibits an alleged similarity to the Cahn-Hilliard equa¬ 
tion studied in the context of liquid-liquid demixing m 
Though both equations have terms of similar orders in 
0(Vp), they are fundamentally different with respect to 
the saturation of droplet or colony growth. The Cahn- 
Hilliard equation exhibits an instability of the homoge¬ 
neous state, which saturates because the effective diffu¬ 
sion constant in front of the LaPlace operator decreases 
to zero. Eq. © also exhibits an instability but it satu¬ 
rates due to a different mechanism as discussed below. 

Colony formation due to pili-mediated interaetion: 
The condition for the instability in Eq. © is a{p) < 0 . 
Its onset marks a critical density, pc = G/pi. Eor 
Po > Pc: fhe homogenous state of density po is unsta¬ 
ble. The instability enhances small density modulations 
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FIG. 2. (color online) (a,b) Critical density pc as a function 
of the non-dimensional parameter G and interaction strength 
a. Each line separates the parameter space, where colonies 
develop or the system remains homogeneous. In (b), three 
values of interaction strength a are displayed: (0.5,0.25,0.1) 
=(solid, dashed, dash-dotted), and (c) depicts three values 
of G: (0.1, 0.2, 0.5)=(solid, dashed, dash-dotted), (c) Illus¬ 
tration of how the instability is balanced: For initial densi¬ 
ties p{t = 0) > Pc the ‘instability flux’ jins drives the emer¬ 
gence of a spatially inhomogeneous density profile. Depend¬ 
ing on the location along the density prohle, the ‘curvature 
flux’ jcu and/or the ‘gradient curvature flux’ jvcu acts against 
the ‘instability flux’ jins and thereby balances the instabil¬ 
ity. (d) Maximal density minus minimal density, pmax — Pmin, 
as a function of time t, where pniax(t) = maxrp(r,t) and 
Pmin(t) — minrp(r, t), for numerical solutions to Ep. (§ with 
and without ‘curvature flux’ jcu- 


around the homogenous density po with a dispersion rela¬ 
tion w{q) = —a{pQ)q‘^ — /^(po)^^. Pc depends on the non- 
dimensional parameter G and the interaction strength a, 
Pc = Gjiaac^. We find that pc decreases for stronger 
attractive interactions, a 1/2, and smaller values of 
G; see Fig. |lla,b). 

The instability is opposed by fluxes related to the spa¬ 
tial curvature of the density field, which can be qual¬ 
itatively understood by splitting the flux, j = jinst + 
jcu + jvcu with dtp = -V • j. jinst = -a{p)Vp denotes 
the ‘instability flux’ which acts for a < 0 like negative 
diffusion thus driving particles to the center of a den¬ 
sity spot [see Fig. |^c) for an illustration]. There the 
instability current is opposed by the ‘curvature flux’, 
jcu = —/d2(V^p)Vp, and the ‘gradient-curvature flux’, 
jvcu = (/^2 + /^3)pV(V^p). Both are directed outwards of 
the density spot since curvature is negative and increases. 


Our findings on the instability and its saturation can 
be scrutinized by numerically solving Eq. Q. A rep¬ 
resentative snapshot of a state at large time-scales is 
shown in Fig. [^b), which appears to be similar to N. 
gonorrhoeae colonies three hours after sedimentation on 
a plastic substrate [Fig. [^a)]. Using parameter values 
consistent with the experimental system we observe mul¬ 
tiple colonies developing quickly for densities above the 
critical value. We checked numerically that for all param¬ 
eter values lying within the ‘colony phase’ of the analytic 
phase diagram [Fig. |^a,b)] give rise to the formation of 
colonies. After the onset of the instability, colonies ex¬ 
ponentially grow with a growth speed that is higher the 
larger the difference of the homogeneous density po to the 
critical density pc; see Supplemental Material [31], S.4. 
Thus, for Po Pc, we observe a colony growth rate 
decreasing to zero; a phenomena reminiscent of ‘critical 
slowing down’ in phase transitions [38]. Subsequent to 
the initial growth, there is regime, where colonies grow 
only very slowly [Fig. [^d), solid red line]. The later 
observation is due to a weak interaction between the 
colonies via some evaporation-condensation mechanism 
qualitatively reminiscent of Ostwald-ripening in liquid- 
liquid phase-separation m Interestingly, at the onset 
of the instability the non-linear ‘curvature flux’ jcu van¬ 
ishes suggesting that it might play an essential role for 
developed colonies at large time-scales. Running the sys¬ 
tem without curvature flux, jcu = 0, we find that the 
subsequent ripening is absent leading to a stable state 
consisting of multiple colonies [see Fig.|^d), dashed line]. 
This implies that interactions between colonies is driven 
by the ‘curvature flux’, while the ‘gradient curvature 
flux’ suffices for the saturation. Based on this insight 
we can analytically estimate the colony size at the time 
when the system crosses to the very slow ripening regime 
[vertical line in Fig. [^d)] by neglecting the curvature 
flux (see Supplemental Material m, S.3 for more de¬ 
tails). For densities p ^ pc^ stationary periodic solu- 
tio ns are suppor ted with the quasi-static colony size of 
7r\/(/d2 + /ds)/;^!. Remarkably, for a = 0.5, this estimate 
suggests a colony size of several pili-lengths, which is con¬ 
sistent with N. gonorrhoeae [Fig. Ba)]. 

Biologieal relevanee: In principle, all parameters en¬ 
tering the kinetic description Eq. 0 can be measured or 
estimated for living colonies forming on a substrate and 
thereby all kinetic coefficients in Eq. 0. In particular, 
for N. gonorrhoeae^ ^ l/rm 01301 and colony forma¬ 
tion is observed for densities of p ~ 0.2. The attachment 
rate to the substrate can be obtained from measurements 
of the single cell diffusion constant, /Cq = ^/(3^pi) ~ 
with D ^ O.b/rm^/s [4] and the cell-cell interac¬ 
tion rate can be roughly estimated from the experimental 
value of the mean next neighbour distance and the mean 
pili-number per cell to Wo/^pi ^ (see Supplemental 
Material [31], S7). Therefore, a typical value for the di¬ 
mensionless parameter for N. gonorrhoeae is G ^ O.lq”^. 
Recently, the attachment probability of pili to a substrate 
has been determined by fitting a model to experimental 
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FIG. 3. (color online) Critical density Pc(<^ad, api)/G as a 
function of (a) adhesive strength a^d (with Upi = 0.5) and 
(b) pili-mediated interaction strength Upi (with a^d = 0.5). 
In both plots, G = 0.1 and each solid line corresponds to 
7 G {0.4, 0.3, 0.2} from top to bottom. Black line corresponds 
to 7 = 0.3. The horizontal black dashed line (plotted only 
for 7 = 0.3 in (a)) marks the dissolution boundary: Below, 
dissolution (red arrow) is possible, above not. Blue shaded 
areas correspond to the “colony-phase”. 


results [9], finding an approximate value of 0.5. We ex¬ 
pect a roughly similar, maybe lower value for 7 since suc¬ 
cessful binding to another cell can be hindered by other 
moving cells. So far an appropriate estimate for the in¬ 
teraction strength a is missing because the synchronous 
visualization of pili and cell movement is not feasible for 
large enough time-scales. Thereby, we consider a as an 
unknown parameter. 

The proposed kinetic description, Eq. Q, can also be 
used to include other attractive interactions such as ad¬ 
hesion. Since cell-cell adhesion constitutes a local inter¬ 
action on the scale of the cell diameter, an appropriate 
weight function is for example a Gaussian of the form 
>Vad(|ri2|/4d) = Wo/(7r£ad)exp (-r?27ad)> where 4d 
denotes the characteristic length scale which is in the or¬ 
der of the cell size. Comparing both interactions (see 
Supplemental Material for details [31], S.5) we find that 
pili allow for a significantly more pronounced affinity for 
colony formation compared to adhesive interactions, i.e. 
colonies already form at smaller initial density of cells. 

Colony dissolution: Many bacteria are known to inter¬ 
act simultaneously by adhesion and pili. It is hypothe¬ 
sized that these bacteria are able to switch off either ad¬ 
hesion or the pili-mediated interaction without affecting 
their ability to move [nnaES]. Now we address the 
question of whether developed colonies can dissolve by 
switching off either one of these interactions. In other 
words, given the phase-diagram of a specific bacteria 
system, we discuss some possible means of leaving the 
“colony-phase” by ^ 0 or Uad —> 0. We now include 
both interactions by adding a term for pili-mediated in¬ 
teractions Cint,pi and a term corresponding to adhesive 
interactions Cint,ad on the right hand side of Eq. (la), i.e. 
Cint = Cint,pi + Cint,ad- In addition to the already intro¬ 
duced different length scales and ^ad^ we also distin¬ 
guish the corresponding interaction strengths, denoted as 


Upi and Uad (values for adhesion and pili-mediated inter¬ 
actions are denoted as c/c,ad and c/c^pi). We rescale co¬ 
ordinates, density and time by the adhesive interaction 
length 4 d (or cell size), i.e. r -> r • 4 d, / ^ //^ad = P 
and t ^ t' ^ad /^0 5 thereby introducing a ratio of these 
length scales, e = ipi/iad- For the case where cells 
interact with both adhesive and pili-mediated interac¬ 
tions, we find the following effective diffusion constant 
(further coefficients see Supplemental Material [31], S 6 ): 
a{p) = G - p [ttadSad C2.ad + ttpiSpi C2,pi 7 e^]. Setting 
this equation equal to zero marks a critical density 
P*c{aad, flpi) = G/ [fladflad C2,ad + OpiSpi C2,pi 7^^] depend- 
ing on the strength of both interactions, Uad and Upi. 

In order to study the impact of both interactions for 
dissolution of colonies we choose the parameters (e, 7 ) 
relevant to N. gonorrhoeae. Eig. [^a) shows p*(aad) as a 
function of Uad for Upi = 0.5 and G = 0 . 1 , while Eig.|^b) 
depicts p*(api) as a function of Upi for Uad = 0.5 and 
G = 0.5, both for several values of 7 . Eor a given 7 , 
there are two qualitatively distinct regimes for the case 
where adhesive interactions are switched off [Eig. ia)]: 
Eor small enough p* below the “dissolution boundary” 
(horizontal dashed line), colonies can dissolve by switch¬ 
ing off the adhesive interaction (uad ^ 0 ) and is indi¬ 
cated by the red arrows. However, above the dissolution 
boundary, colonies cannot dissolve. Interestingly, choos¬ 
ing the parameters relevant to N. gonorrhoeae gives a 
rather small density regime, where colonies can dissolve, 
rendering the dissolution scenario through switching off 
adhesion as a non-robust mechanism. This is in stark 
contrast to the scenario of switching off pili-mediated in¬ 
teractions [Eig. [^b)]: Eor a given 7 , dissolution is pos¬ 
sible for all experimental densities in the “colony-phase” 
by lowering the pili-interaction strength, Upi ^ 0. These 
findings suggest that switching off pili-mediated interac¬ 
tions is a more robust mechanism for the dissolution of 
bacterial colonies than switching off adhesion. 

To summarize, the formation of living colonies is in¬ 
vestigated using a hydrodynamic equation derived from 
a kinetic description, where most of the parameters can 
be estimated from experimental data for N. gonorrhoeae 
bacteria. Our results demonstrate that kinetic theory can 
be successfully used to describe complex far from equi¬ 
librium systems such as formation and dissolution of liv¬ 
ing bacterial colonies. Applications of this theory could 
pave the way for the physical quantification of the initial 
stages of biofilm formation. Though biological reasons 
for colony formation are specific to each system there 
are qualitative similarities m- Colonies form due to 
encounters with nearby individuals giving rise to struc¬ 
tures of a characteristic size determined by the intra¬ 
species interactions and the environment. These similar¬ 
ities suggest that our kinetic description might be applied 
to other colony-forming systems while the kinetic coeffi¬ 
cients in the resulting hydrodynamic equation may dif¬ 
fer for each system. Eurther open questions concern the 
role of cell division and stochastic fluctuations in living 
colonies HOI. 
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SI. Coordinate change due to isotropy of 
interaction kernel 

Using the collision rule (see Eq. ( 2 ), main text; and 
Fig. SI), 

(ri,r2) ^ ( 7 ^l, 7 ^ 2 ) = (ri - a ri2, r2 + a ri2) , (SI) 
the gain term can be written as: 

= y dri J dr2 W(ri,r2)/(ri,t)/(r2,t) 

^ ((^i - a ■ ri2) - r) + (5 ((r2 + a ■ ri2) “ ^ 

(S2) 


Similar manipulations can be performed for the loss term: 
c~ = j dri J dr 2 W(|ri2|)/(ri,i)/(r2,t)(5(r2 - r) 

= J dri2 J dRi2W(|ri2|)5(Ri2-(ri2/2 + r)) 

X /(Iti2 + ri2/2, t)/(Ri2 — r 12 12 ,t) 

= ydri2 W(|ri2|)/(r+ ri2) 

= f{r,t) /'dri2 W(|ri2|)/(r - ri2). 

(S6) 

Please note that non-local integrands above resemble a 
phenomenological description for the assembly of active 
bundles [ 1 - 3 ]. 


The equation above can be rewritten in terms of relative 
coordinates ri2 = ri — r2 and center-of-mass coordinates 
Ri 2 = (i*i +1*2)72, i.e. ri = Ri 2 + ri2/2 and r2 = R12 - 
T12/2 (see Fig. SI): 


dri2 y dRi2W(|ri2|) 

X /(R-12 + ri2/2, t)f (Ri 2 - ri2/2, t) 


S ( Ri 2 + ri2 ( - 


+ <5 ( Ri 2 + ri2 


y+a 


(S 3 ) 


S2. Details of coarse-graining and truncation 


Truncating of the non-local distribution function /(r± 
ari2) (cf. Eq. (S 4 ) and Eq. (S6)) at the fourth order leads 
to 


/(r ± ari2) = /(r) ± a (ri2 • V) /(r) + ^ ■ Vf /(r) 


± Y (^^12 • V) /(r) + — (ri2 ■ V) /(r) + O [(ri2 • V) 


5 

(S 7 ) 


Due to the isotropy of the collision kernel the integration 
over the center-of-mass coordinates can be performed, 
finding 




= dri2 W(|ri2|) 


/(r + ari2,t)/(r-ari2,i) 


+ f{r-ari2,t)f{r + ari2,t) , 


(S 4 ) 


where a = 1 — a. The equation above can be further 
simplified for example for the case a = ^: 


C+= 7dri2 W(|ri2|)/(r + ri2/2,i)/(r-ri2/2,i). 

(S 5 ) 



FIG. SI. Illustration of collision setup and the collision 
rule mimicking attractive interactions, (a) The two collision 
partners have the spatial coordinates ri and r 2 , defining a 
relative distance ri 2 = ri — r 2 and a center-of-mass coordinate 
Ri 2 = (ri -h 1 * 2 ) 72 . (b) The collision given in Eq. (2) [main 
text], which maps the pre-collision coordinates (ri,r 2 ) to the 
post-collision coordinates (IZi (ri,r 2 ) , 7^2 (ri,r 2 )) [indicated 
by green arrows]. For the illustration, a = 0.25. 
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where we omitted the time dependence for reasons of 
brevity. Using the truncation above and neglecting 


all terms O 


(1*12 • V) amounts to an explicit coarse- 


graining of the system’s dynamics to length scales beyond 
the characteristic length scale of the interaction. There¬ 
fore, we will refer to the resulting equation obtained at 
the end of this section as hydrodynamic equation. 


1. Single cell motility term: 


The term modeling the single cell motility across the 
substate, Cmot? can be coarse-grained as follows. Defining 
the pili-mediated displacement as b = (bx^by) = r' — r, 
and the transition rate /Cr^r' = A^(b; r): 

Cmot(r,t) = J dr' [ICr>^rf{r',t) - ICr^r' 

= y db [/C(-b; r + b)/(r', t) - JC{h; r)/(r, t)] 

= j db[/C(b;r-b)/(r-b,i)-/C(b;r)/(r,t)] 

(S8) 


Expanding the non-local integrand with respect to the 
spatial coordinates and keeping only the highest non¬ 
vanishing order leads to: 


Cmot(r, ^) = y (b • V)^ [/C(b; r)/(r, t)] 

= y dbi/C(b) {bxdx + bydyf f{r, t) 
= y dbi/C(b) [bldl + bpl) f{r,t) 

= JdhlC(\h\)bl {dl + d^y)fir,t) 

= DVV(r,i), 


(S9) 


where we assumed that the transition rate does not de¬ 
pend on the spatial coordinates, /C(b;r) = /C(b), and 
that it is an even [4] and isotropic function, /C(b) = 
/C(|b|). Moreover, using /C(|b|) = [-\h\ypi) 

pi 

one obtains the diffusion constant for pili-mediated motil¬ 
ity: 


D = pi, (SIO) 

where C2,pi given in Table SI and /Co denotes the attach¬ 
ment rate of pili to the substrate. Restriction to the 
second order is validated by the Pawula theorem [5]. 


2. Interaction term: 

The interaction term Cint can be split in a gain term, 
C+, and a loss term, C~. 


Gain: Neglecting all terms above the fourth order, we 
find for the loss term: 


C = /(r) co/(r) -f 

where the coefficients are given as 

Cfe = y y dri2,a=dri2,j, r'^2,x bVad/pi(ri2) 

f^/pi-Cfc if fc even, 

I 0 if k odd. 


(Sll) 


(S12) 


with ri2 = |ri2| = y ^i2,cc + ^12,?/ are dimen¬ 

sionless numbers given in Table SI for the adhesive(ad) 
and pih-mediated(pi) interaction. Note that odd pow¬ 
ers, thereby also mixed gradients such as dxdyf, vanish 
since the integral is then an asymmetric function with 
respect to the integration over e.g. dri 2 ^x or dri2,cc, re¬ 
spectively [6]. 

Loss: Following similar lines for the source term C+, 
and neglecting all contributions above the fourth order 
in the spatial derivatives, we find only six non-zero con¬ 
tributions: 

Cint = Co + ^[a^ + d^]c2 /(r)VV(r) - aac2 |V/(r )|2 

+ + d\4 + ^a^a‘^C4 

-\[aa^+da^]c4 [V/(r)] • Vy(r). 

(S13) 


Gain and loss: Combining gain and loss term, Cint = 
C+ — C“, and plugging it in Eq. (1) [main text], leads to 
the final hydrodynamic equation: 

dtfir, t) = -DV^/(r, f)] + - l]c2 /(r, f)V^/(r, t) 

+ ^[a^ + a^ - 1]C4/(r,f)V'‘/(r, t) 

-adC2 |V/(r,f)|^ -h ^a‘^d‘^C4 [V^/(r,f)]^ 

- \ [ad^ +da^]c4 [V/(r,f)] ■V^f{r,t). 

(S14) 


Note that in Eq. (S14) the zeroth order term cancels be¬ 
cause of particle conservation. 

As last step, we use the scaling of the kinetic coeffi¬ 
cients, Ck = Woi^Ck [see Eq. (S12)], with i G {^ad,'^pi} 
and the numerical values given in Table SI, and write 
Eq. (S14) in a dimensionless form. To this end, we rescale 
the coordinates and the density by means of the interac¬ 
tion length i.e. r ^ rT and f ^ f /P = p. This implies 
a rescaling of the time-scale given by t ^ 
note that 7 = 1 for adhesive interactions. Along these 
lines we find the following dimensionless parameter 


7W0’ 


(S15) 



3 


k 

0 12 3 4 

ek,a,d 

10 i 0 f 

Cfe,pi 

1 0 3 0 45 


TABLE SI. The numerical numbers Ck (as defined in 
Eq. (SI2)) corresponding to adhesion and pili-mediated inter¬ 
actions (refer to main text for the respective collision kernels). 


where D denotes the single cell diffusion constant given 
by Eq. (SIO). 

Using the aforementioned rescalings leads to the fol¬ 
lowing dimensionless equation: 

dtp{r, t) = «(/>) V V(i’, t) + k{p) V‘/>(r, t) - jSi | V/)(r, t)f 


+ P 2 [VV(r, i)] ^ -/33 [Vp{r,t)]-V^p{r,t) , 



(S16) 

where the kinetic coefficients are defined as: 


a{p) = G- 01 p(r,t), 

(S17) 

k{p) = - (02 + 03) /3(r,t), 

(S18) 

01 = aac2, 

(S19) 

to 

II 

h- 

to 

to 

(S 20 ) 

,^3 = i (aa^ + aa^) £ 4 , 

(S 21 ) 


and with numerical values given in Table SI and a = 
1 — a. It is worth mentioning that Eq. (S16) conserves 
density, i.e. it can be written as dtp = —V • j, with the 
flux j given by [7]: 

j = -a{p)Vp - [k{p)V^ + 02 (VV)] Vp. (S 22 ) 


S3. Quasi-stationary colony size 

An approximative expression for the colony size when 
the system transits from fast initial growth to the slow 
ripening regime (see main text) can be obtained by 
neglecting the ‘curvature flux’, jcu = P 2 (V^p)V/) in 
Eq. (S 22 ). Note that the ‘curvature flux’ also vanishes af¬ 
ter linearization of Eq. (S 22 ) around po with p = pQ^dp. 
Assuming quasi-static conditions, j = 0 , leads to 

0 = [a{po) + k(po)V2] V5p. (S23) 

Within this quasi-static approximation there are two sta¬ 
tionary states: The homogenous density field p = po with 
Vp = 0 and an inhomogeneous state (Vp 7 ^ 0) that sup¬ 
ports periodic solutions suggesting the coexistence of sev¬ 
eral colonies. Since our simulations indicate that droplets 
have a very similar size when the system crosses from the 
fast initial growth to the slow ripening phase (see video 
material), let us approximate the droplet size distribu¬ 
tion to be infinitely narrow and extract a single length 


scale, referred to as quasi-stationary colony size Writ¬ 
ing V ^ iq and ^ = 7r|q|“^, one finds for Vp 7 ^ 0 : 

O:{po) 

For large densities, one gets: 

C(po 00) ~ 

which gives approximately 5 pili-length for a = 0.5; a 
value that is consistent with N. gonorrhoeae colonies 
forming on a substrate [see Fig. 1(a), main text]. 


S4. Numerical solution 

Eq. (SI 6 ) was solved numerically using a 8 th (embed¬ 
ded 9th) order adaptive Runge-Kutta method. Specif¬ 
ically, we employed an established package called 
“XMDS 2 ” (see [ 8 ]) developed by Graham R. Dennis, 
Joseph J. Hope and Mattias T. Johnsson (Comput. Phys. 
Commun. 184, 201-208 (2013)). 

Fig. S 2 depicts numerical solutions of Eq. (S16) for 
various density values in order to illustrate that the 
equation—originating from a gradient expansion and a 
truncation at the fourth order—is well-behaved in case 
of an instability. For a discussion of the results, please 
refer to the figure caption of Fig. S 2 . 


S5. Phase diagrams: Pili-mediated interaction 
versus adhesion 

Since the pili length distribution exhibits a more pro¬ 
nounced tail than the localised Gaussian distribution 
(characterised by larger C 2 ) and also a larger characteris¬ 
tic length as for N. gonorrhoeae > ^ad)? pili allow for 
a significantly more pronounced affinity for colony for¬ 
mation compared to adhesive interactions [Fig. S3(a,b)]. 


S6. Kinetic coefficients for adhesion and 
pili-mediated interactions 


For adhesive and pili-mediated interactions, solely the 
effective diffusion constant a{p) has been given in the 
main text. Reading the coefficients Eq. (S19)-(S21) as 
function of the interaction strength a, i.e. Pi = Pi{a)^ the 
remaining coefficient is listed below: 


K,{p) 


{02{aa.d) + 0i{aa.d)) 


+ {02{api) +/33(api)) 7 


(S26) 


The corresponding length scales, interaction strength 
and numerical coefficients for adhesion and pili- 
mediated interactions are denoted as: ^pi, c/c,ad and 
4d,aad,4.ad, respectively, and e = 4i/4d- The rescal- 
ings used are described in the main text. 
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FIG. S2. (color online) (a,b) Maximal density minus min¬ 
imal density, pmax — Pmin, as a function of time t, where 
Pmax(t) = maxrp(r,t) and pmin(t) = minrp(r,t). In (a) 
we divide through the mean density, po = J drp(r,t), 
where L is the system size (here we used L — 20). Note 
that for all numerical runs, p is conserved. For the numer¬ 
ical results (a-d) we used a = 0.5, G = 0.1 and consider 
pili-mediated interactions (c 2 = 3/4), implying a critical den¬ 
sity of Pc = 2/15 = 0.133 (see main text for the analytic 
expression of the critical density), (a) Consistently, realiza¬ 
tions with p > Pc exhibit an instability, i.e. pmax — Pmin grows 
as a function of time, whereas for po < Pc, weak initial spa¬ 
tial perturbations around po decay exponentially. The growth 
roughly follows an exponential with a speed strongly depen¬ 
dent on the difference of po to the critical density: The larger 
this difference, the faster the initial growth speed. Please 
note that for densities slightly above the critical threshold 
(e.g. Po = 0.2), the instability grows too slowly to capture 
the long time behavior. However, for large enough densities 
(e.g. Po = {0.4, 0.6,1.0}), Pmax — Pmin clearly indicates a sat¬ 
uration. (b) The behavior at very large time-scales is very 
hard to capture numerically: Interestingly, if the system has 
developed to a single colony, pmax — Pmin becomes ffat (see 
Po = 0.4 and snapshot (c)). If there are two or more colonies 
in the system (e.g. po = {1-0} and snapshot (d)), pmax — Pmin 
still changes as a function of time, though very weakly. The 
details of this coarsening mechanism close to saturation will 
be studied elsewhere. 


S7. Estimates of parameters for N. gonorrhoeae 

1 . Estimate for density po 

We considered 20 experimental realizations of N. gon¬ 
orrhoeae bacteria forming colony in a substrate. Af¬ 
ter sedimentation to the plastic surface we calculated 
from the corresponding binary images the overall area 




FIG. S3, (color online) (a,b) Gritical density pc as a function 
of the non-dimensional parameter G and interaction strength 
a for pili-mediated interactions (blue) and adhesion (red). 
In (a), three values of interaction strength a are displayed: 
(0.5, 0.25,0.1) = (solid,dashed,dash-dotted), ipi ^ 2 ■ and 
in (b), G = 0.1. 


fraction covered by bacteria cells, finding (j) ^ 0.1. 
The corresponding dimensionless density is then: po = 
(Tri^ceii) ~ 0.125 with = Ijim and Rcew = 0.5pm. 
Since typically some small three dimensional colonies 
have already formed during sedimentation process we ex¬ 
pect that the determined value represents a slight under¬ 
estimation, thereby we use po = 0-2 in the manuscript. 


2. Estimation of Wo f ipi 


is a measure for the rate of cell-cell encounters 
occurring in an area of ^p-. A direct measurement of this 
quantity requires the sampling of the cell trajectories on 
the time-scale of the cell-cell encounters which is intricate 
because bacteria cells get harmed in case of too frequent 
light exposure. Therefore, we have to content with a 
rough estimate. If we assume that each pilus per cell 
acts independently, we can first estimate the rate of a 
cell-cell interaction for a single pilus. For intermediate 
and large cell densities with respect to the intersection 
scales, it is expected that this rate roughly scales with 
the number of pili per cell. 

The time r between two interactions using a single 
pilus should be roughly given by the time to diffuse 
the distance to the next-neighbouring cell ^nn, ^ ^ 
i‘^^/{4:D), where D = 0.5pm^/s. However, this distance 
is reduced by the cell-diameter 2Rce\\ = and two 

times the typical pili-length 2ip[ = 2pm [9], leading to: 

r ^ (^nn — 2 i?ceii — 2^pi)^ /(4^) ~ 25, where we deter¬ 
mined ^nn ~ 5pm from the binary images directly after 
sedimentation. 

TEM-images have revealed a mean pili number in or¬ 
der of A" ^ 10 [10]. For intermediate and large cell 
densities, the overall-interaction rate for for N. gonor¬ 
rhoeae at surface coverage of (j) ^ 0.1 is approximately 
Wo/^pi N . (25)-i 55 -1. Thereby, G ^ 0.1. 
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S8. Video Material 


Using the procedure mentioned in Section S4 we solved 
the hydrodynamic equation for various parameters sets. 


For a representative parameter set (G = 0.1; Uad =0.0; 
ttpi = 0.5; po = 0.6), a video is appended (T denotes 
the dimensionless time; see main text for corresponding 
rescaling): 

C olony-f ormationjvia.kinetic-description.mpA. 
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